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Spin-orbit coupled bosons can exhibit rich equilibrium phases at low temperature and in the 
presence of particle-particle interactions. In the case with a ID synthetic spin-orbit interaction, 
it has been observed that the ground state of a Bose gas can be a normal phase, stripe phase, or 
magnetized phase in different experimentally controllable parameter regimes. The magnetized states 
are doubly degenerate and consist of a many-particle two-state system. In this work, we investigate 
the nonequilibrium quantum dynamics by switching on an external perturbation to induce resonant 
couplings between the magnetized phases, and predict the novel quantum spin dynamics which 
cannot be obtained in the single-particle systems. In particular, due to particle-particle interactions, 
the transition of the Bose condensate from one magnetized phase to the other is forbidden when 
the strength of external perturbation is less than a critical value, and a full transition can occur 
only when the perturbation exceeds such critical strength. This phenomenon manifests itself a 
quantum dynamical phase transition, with the critical point behavior being exactly solvable. From 
the numerical simulations and exact analytic studies we show that the predicted many-body effects 
can be well observed with the current experiments. 

PACS numbers: 


Introduction 

In quantum mechanics, a resonant quantum oscillation 
can occur between two single-particle states when the 
system is subject to external perturbation which drives 
resonant couplings between the states [T]. The full tran¬ 
sition from one quantum state to another can occur after 
one-half period evolution, no matter how weak the res¬ 
onant coupling is. Such two-state mechanism has been 
at the heart of understanding fundamental quantum dy¬ 
namics and developing applications in atomic, molecular, 
optical, and condensed matter systems. A Bose-Einstein 
condensate (BEC) of many atoms can be described with 
a global wave function, which resembles a single-particle 
quantum state [2]. In the absence of particle-particle in¬ 
teractions, the quantum dynamics of a BEC are similar 
as those in a single-particle system. 

Recently, a breakthrough progress for cold atoms is the 
experimental realization of synthetic spin-orbit (SO) cou¬ 
pling for (pseudo)spin-l/2 bosons [3] and fermions [HIS], 
which has generated exciting pursuit of exotic physics 
including the novel SO effects EHIH] and topological 
phases [T3H37] , and open new possibilities to probe novel 
quantum spin dynamics. So far the realized SO interac¬ 
tion is of the one-dimensional (ID) form, and is generated 
through two-photon Raman process which couples the 
two spin states and transfers momentum between them 
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simultaneously [7l [TO]. The SO coupling brings a rich 
structure in the single particle spectra of which the low¬ 
est energy subband can be in the single-well or double¬ 
well form by tuning the Raman coupling strength m- 
For bosons, the double-well dispersion relation is partic¬ 
ularly interesting, for it has two degenerate minima with 
opposite finite momenta in the single-particle spectrum. 
As a consequence of the competition between the kinetic 
and interaction energies, at low temperature the bosons 
may condense in a fixed finite momentum state (magne¬ 
tized phase), which is of two-fold degeneracy, or in the 
superposition of two degenerate states with opposite fi¬ 
nite momenta (stripe phase) [28H3T| . These novel phases 
have been recently observed in the experiment with ®^Rb 
atoms [33]. Rather than being single-particle states, the 
two magnetized phases with opposite finite momenta con¬ 
sist of a degenerate many-particle two-state system. This 
brings about an interesting question that what nontrivial 
dynamics can be expected by applying an external per¬ 
turbation to couple resonantly such two many-particle 
ground states? 

Results. -In this work, we investigate the transitions 
between different many-particle states of a SO coupled 
BEC induced by external perturbation, and uncover 
novel quantum spin dynamics which cannot be obtained 
from single-particle systems. The main results are out¬ 
lined below. Starting from the magnetized phase with 
the momentum kj. = km, an external perturbation is 
switched on to drive transition from the initial phase to 
the one at fc = —km- (i) A critical value Vp^cnt, which 
depends on the particle-particle interactions, is predicted 
for the perturbation. When the strength of perturbation 
is below the critical value, the maximum transition ra¬ 
tio, given by l/S^axL i® always less than 1/2. (ii) By 
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increasing the perturbation strength to be above t^,crit, 
we find that the maximum transition ratio sudden jumps 
from |/3maxl < 1/2 to I^S^axI = 1; rendering a novel quan¬ 
tum dynamical phase transition |33H35j . (iii) The time of 
transition between the two ground phases diverges log¬ 
arithmically as the perturbation approaches the critical 
point, while in the case with the perturbation far above 
the critical magnitude, the transition becomes similar to 
the quantum dynamics of a single-particle two-state sys¬ 
tem. These predictions can be observed in the experi¬ 
ment by measuring the time evolution of spin polariza¬ 
tion of the condensate. 

Magnetized Phases and Induced Couplings 

The system we consider is a spin-1/2 BEC of ®^Rb 
atoms with ID synthetic SO interaction generated by 
two-photon Raman coupling I1I3S1- In the presence of 
the particle-particle interactions, the Hamiltonian of the 
system is given by (let h = m = 1 for simplicity) 

Hbe = + 

where the single-particle term Hq includes ID SO cou¬ 
pling along X direction and is the interacting term 

Ho = ^-f 

Hint = y (1) 

Here the field operator (V'J (i")) annihilates (creates) 
a particle at position r and in the spin state s =t,4., 
the two-photon Raman coupling strength O serves as 
an effective Zeeman field along x axis, and the SO co¬ 
efficient fco is determined by the wave vector of Raman 
lasers. The interacting coefficients gs = and 

9a = j where gss' represents the interacting coefficient 
between atomic species s and s'. We have neglected the 
difference between g^^ and g^, since it can be absorbed 
as an effective two-photon detuning in the Raman cou¬ 
pling and be compensated by properly tuning the fre¬ 
quencies of Raman lasers in the experiment m- Diag¬ 
onalizing the single-particle Hamiltonian Hq yields two 

subbands with energies ± -p 1^. The 

lower subband energy E^^^ has two degenerate minima 
at ifciuin = ±fco(l ~ D^/4fcg)^/^ when |r2| < 2fcg, while 
it has a single minimum at kmin = 0 if |D| > 2kQ [TU] . 
The former band structure leads to nontrivial phases for 
bosons even at zero temperature. 

The Hamiltonian Hbe lias three distinct ground phases 
due to the competition between the kinetic and inter¬ 
action energies which can be controlled in the exper¬ 
iment by manipulating Q and the condensate density 


n = -|- n^, where represents the atomic density of 

the spin component s. In general, the condensate wave 
function of these ground phases take the following form 

$BE(r) = Vn[auR{e)E'^-'^ + /3ML(0)e-*'=™"], (2) 

where ur = (cos 0,— sin 0)"^ and ur = {sin 9, — cos 9)'^ 
with tan20 = H/(2fcofcm) are two single-particle eigen¬ 
states of the lower subband of Hq with momenta kx = 
±km (generically different from ±fc„iin) [see Fig. 0(a,b)], 
a, P and km are variational parameters and should be de¬ 
termined by minimizing the total energy of the system. 
The three different phases are separated by two criti¬ 
cal Raman strengths Hci,c 2 , which depend on condensate 
density and satisfy Hd > m- In the large Raman 
coupling regime, namely, when H > Hd, the total energy 
minimizes at km = 0, giving the trivial condensate phase. 
In the small Raman coupling regime with D < 12^2, the 
ground phase is a superposition of the states with finite 
momenta kx = Ekm and jaj = |/3| = resulting 

in the stripe phase. Finally, for the intermediate regime 
with r 2 c 2 < H < Hci, the system has two degenerate 
phases at finite momenta kx = km and kx = —km, with 
a = l,/3 = 0 and a = 0,/3 = 1, respectively, called the 
magnetized states. These novel phases have been recently 
observed in the experiment with ®^Rb at temperature of 
a few tens of nK [5^ . 

The magnetized phases break both the U{1) and Z 2 
mirror symmetries of the Hamiltonian Hbe, with the mir¬ 
ror symmetry defined by A4x = ic^x and connecting the 
two magnetized phases at kx = Ekm- The main purpose 
of this work is to study the quantum spin dynamics for 
this two state system by introducing an external pertur¬ 
bation to couple resonantly the two magnetized phases 
[see Fig. (c)]. Note that the magnetized phases dis¬ 
tinct from each other by a finite momentum 5k = 2km- 
To induce quantum transitions between such two phases, 
we switch on a perturbation at the time f = 0 in the 
following form 

Jo, t<0 

I4x(r,t) “ I J (prcos^ kmX + V’jV’;) , t > 0. 

where Vq represents the perturbation strength. The 
above perturbation is simply an additional lattice poten¬ 
tial which can be readily generated by a standing wave 
light with wave vector kR = km- We note that in general 
Vex can induce couplings between any two states with 
momentum difference E2km- Thus it may drive transi¬ 
tions from the magnetized phases to other states with e.g. 
the momenta kx = ±3fcm- These states, however, have 
large energy mismatch compared with the energy of the 
ground phases. Taking into account such transitions does 
not affect the main results of this work, as we shall discuss 
later. Below for convenience we first neglect these addi¬ 
tional couplings. On the other hand, from Fig.[T][a,b) we 
read that for magnetized phases km — fcmin — which 
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FIG. 1: The condensate momentum and induced cou¬ 
pling between magnetized phases, a. The ratio between 
the condensate momentum km and fco as a function of Ra¬ 
man coupling strength Q with different interaction energies 
gall (in the unit of fep)- b. The ratio between the condensate 
momentum km and fcmin- As the Raman coupling Q. increases, 
the ratios km/ko and km/kmm increase for the stripe phase 
regime, decrease for the magnetized phase regime, and is zero 
for the normal regime, respectively, c. Starting with a mag¬ 
netized phase, an additional external potential switched on at 
t = 0 can drive resonant couplings between the two degener¬ 
ate magnetized phases. 


can facilitate the parameter choice in the real experiment 
by setting km = ^min or kji = Uq. The small mismatch of 
the wave vector only has a tiny quantitative correction 
to the main results, and will be examined later. 

Assume that the system is initialized in the mag¬ 
netized phase with k^ = km- The perturbation I4x 
drives the system to evolve from one magnetized phase 
to another. Keeping only the transitions between the 
two ground states, we can describe the time-dependent 
condensate wave function as |fE'BE(r,t)) = -b 

f3* \vac), where N is the number of atoms of the 

BEG, \vac) represents the vacuum state, and the field op¬ 
erator ijjfi = (cos^V't ~ and V'l = (sin0i/)^ — 

coswith the initial condition a(0) = 1 and 
/3(0) = 0. Denoting by = a{t)ipR + (3{t)ipL, the dy¬ 
namics of the condensate can be derived according to the 
following equation: 

= m),Ho + H-mt + Kx]. (4) 

The right hand side of the above formula contains three 


commutators, corresponding to the kinetic energy term, 
the interacting energy term, and the coupling induced by 
external perturbation between two magnetized phases, 
respectively. It will be shown that due to the particle- 
particle interactions the quantum spin dynamics for the 
many-particle ground states are fundamentally different 
from those in the single-particle systems. 


Quantum Oscillation of Magnetized Phases 

Equation of Motion. -We first derive explicitly the 
equation of motion for the condensate by keeping only the 
transitions between the two magnetized states. This can 
be performed by calculating the commutators in Eq. (|^ 
and projecting the results onto the ground states gener¬ 
ated by iliR and tpL- The resulted equation in the mean- 
field approximation governs the quantum spin dynamics 
for the magnetized phases. The relevant commutators 
can be verified by (details can be found in the Supple¬ 
mentary Material [HH] ! 

[i’R/L,Ho] = 

['4’R/L,Ve^] = iVol2)lpRiL + VptpR/R, 

[^R/L,Hmi] =[G,±Emi\a\^-m^)]^I^R/L 

where Vp = ^ cos 9 sin 9 represents the coupling strength 
induced by external perturbation, Gi = j{gs+ga)j G 2 = 
jigs ~ go), Em = G 2 cos^ 29, and Es = 2Gi cos^ 0sin^ 9. 
Erom the last commutator we can see that when a/3 0, 

the interacting Hamiltonian also contributes to the cou¬ 
pling between ipR and ipR. This is because in this case 
the condensate wave function is a superposition of the 
ground states at kx = Ekm, which breaks translational 
symmetry and leads to a periodic density prohle with 
spatial period Sx = ir/km in the condensate. Accord¬ 
ingly, the interacting Hamiltonian E[mt also breaks trans¬ 
lational symmetry and can drive transitions between the 
two ground phases by transferring momentum 2km be¬ 
tween them. 

With the above results we can now write down the 
equation of motion in the space of ipR and 'ipR- 

where the effective two-state nonlinear Hamiltonian is 
given by 


'kLefi — E^^^ + ^ + Gi + Emi\a\'^ — \l3\‘^)sz 

+2Ea [3?(a/3*)S, - ^iap*)sy] + VpSx- (7) 

Here Sx,y,z are Pauli matrices acting on the pseudospin 
space spanned by the two magnetized states. It can be 





















4 


seen that Em and Eg in the above formula characterize 
the interaction energies for the magnetized phase and 
stripe phase, respectively. Actually, it is easy to verify 
that without external perturbation, the magnetized 
phase, say with a = 1 and /3 = 0, is an eigenstate of HeS 
with energy E = Eq + Em, where Eq = E~^ + ^ + Gi, 
and the energy of a stripe state with |a| = |/3| = l/-\/2 
is E = Eq + Eg. The BEG is initialized in one of the 
magnetized phases if Em < Eg, while in the opposite 
case the BEG will be initialized in a stripe phase. 

Pseudospin Evolution.- We numerically solve the 
nonlinear equation of motion with the initial condition 
a(0) = 1 and /3(0) = 0. The real spin evolution can be 
easily obtained from its relation to the pseudospin dy¬ 
namics (o'z(t)) = (sj;(<)) COS 20, with {sz{t)) = |ap —|/?p, 
= sin20/2, and (tfy) = 0. Eig. shows the 
pseudospin evolution, from which one can determine 
that time-dependent probability, given by |/3(t)P = [1 — 
(sz(t))]/2, of the magnetized state at kx = —km with 
different perturbation amplitudes and under experimen¬ 
tally achievable parameter regimes. Erom the numerical 
results, we find the following novel features. 

Eirst of all, a minimum value of the perturbation 
l^,crit, which relates to the critical amplitude of 14x by 
ho.crit = EVp^cmt/ sin 20 and depends on the condensate 
interaction energy ggn, is required to have a full tran¬ 
sition from one magnetized phase to another. When 
the perturbation strength is below this critical value, the 
pseudospin evolves only within the upper half spherical 
surface and the maximum probability of the magnetized 
phase at kx = —km is less than 1/2, i.e. |/3maxl < 1/2 
[Eig. I (a,b)]. On the other hand, once the the perturba¬ 
tion exceeds the critical value, the pseudospin can pass 
through the equatorial {x — y) plane, where the conden¬ 
sate is in a stripe state, and the maximum transition ratio 
jumps from I/S^axl < 1/2 to |/3maxl = 1> thus a full tran¬ 
sition can occur [Fig. (c,d)]. From Fig. (b,c) we can 
see a sudden change in the pseudospin evolution (i(t)) 
from Vp < Ip.crit to 1^ > Ip.crit, showing the interesting 
critical behavior of dynamics for the transitions between 
the two many-particle ground states. 

It can be seen that the period of pseudospin evolution 
Tr exhibits an increasing function of the perturbation 
when Vp < Ip.crit (or Vq < Vb^crit), and a decreasing 
function of the perturbation when > V/_crit- Moreover, 
from Fig. Kb, c) one can read that the pseudospin evolves 
faster around (sz) ~ ±1 than it does around (Sz) ~ 0. 

The existence of the critical perturbation reveals a 
novel quantum dynamical phase transition for the many- 
particle states driven by the applied perturbation [31- 
|35l [38] , which cannot be obtained from a single-particle 
two-state system. This phenomenon reflects the nontriv¬ 
ial interactions effects on the transitions between many- 
particle states, and can be qualitatively understood in the 
following way. When a quantum system evolves from one 


state to another, it must undergo the intermediate pro¬ 
cess which is a superposition of the two quantum states. 
For single-particle system with resonant couplings, the 
superposition of the two states has exactly the same en¬ 
ergy as any of the two quantum states. Thus a resonant 
oscillation is guaranteed in the whole process, and a full 
transition can always happen as long as the perturba¬ 
tion is nonzero. However, for the present many-particle 
interacting system, the transition from one magnetized 
phase to the other must undergo the intermediate stripe 
phase which has the energy different from (greater than) 
that of the magnetized phases. As a result, the transition 
between the two ground phases is “detuned” from the in¬ 
termediate superposition state, and then a full transition 
is forbidden if the coupling is not strong enough. 


Effective Theory for Pseudospin Dynamics 


Now we turn to the analytical studies on the quan¬ 
tum spin dynamics. We shall develop an effective theory 


based on the equation of motion (S51, with which we can 


reach a full understanding of the results obtained in the 
previous section. It is convenient to recast the equation 
of motion into a set of new formulas for the real variables 
introduced by: = \a\’^,i 2 = Q'(q:*/3),^3 = ^{a*l3). It 

follows that d^i/dt = 2 Vp^ 2 , d^ 2 /dt = (1 — 2 ^i)(I^-|-A^ 3 ), 
and d^s/dt = —2A(1 — 2^i)^2, with A = Eg — Em- Note 
that the initial condition is ^i(O) = 1 and ^ 2 ( 0 ) = ^3(0) = 
0, from which we further show that these equations can 
be organized and simplified into the following dynamical 
equations (see Supplementary Material for details [55]1: 


d2(sz(t)) , dV;ff(sz)) 


dt'^ 

{Sy{t)) = 


d{s,,) 

dt{sz{t)) 


= 0 , 


{Sx{t)) = — [ 


2Vp 

1 rdt{Sy{t)) 


( 8 ) 


A^2{sz{t)) 

where Vbff = {2Vp - A^){sz)'^ E (A^/2)(sz)'‘ is a 
type function of (sz) and can be in the double-well (for 
A > ^/2Vp) or single-well (for A < V2Vp) form. The 
above equations provide a clear interpretation for the 
quantum spin dynamics of the present SO coupled in¬ 
teracting bosons. Namely, the spin evolution (Sz(t)) can 
be effectively treated as a classical particle which moves 
in the effective potential I4ff((sz)), with the initial posi¬ 
tion ( 32 ( 0 )) = 1 and initial velocity d(sz)/dt|t=o = 0 [see 
Fig.[^(a,b)]. From the formula Q, we can exactly study 
the quantum transitions between the magnetized phases, 
and obtain the following important results. 

First, the critical external perturbation can be read 
from the condition 14ff(l) = 14ff(0), which follows that 


l/,,crit — 2 — 2^^^ ^m)- 


(9) 
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FIG. 2: Pseudospin spin evolution on the Bloch spherical surface. The arrows depict the direction of pseudospin {s{t)} 
during a single evolution period which starts from the north point, with the color quantifying the evolution time. a. The 
external perturbation is below the critical value, with Vo = 0.6Vb,crit, where Vb.crit = 41p,crit/ sin 20. The evolution path is 
nearly a circle for the small perturbation regime, b. For the external perturbation below but close to the critical value, with 
Vo = 0.999Vo,crit, the evolution path exhibits a “heart” shape. At the southmost point the pseudospin vector points closely 
to —X direction, and corresponds to the maximum probability value |/3maxl ^ 1/2 of the magnetized phase at k = —km- c. 
The pseudospin evolution with Vo = l.OOlVb.crit. A full transition from the initial magnetized phase to the one at k = —km 
occurs in this case. The evolution path exhibits an “8” shape, and pass through the equatorial {x — y) plane through the point 
close to —X axis. d. The perturbation Vo = 1.4Vo,crit. Further increasing the perturbation broadens the evolution path. Other 
parameters are taken as gsU = l.Ofcg and Q = O.Sfcg in the numerical simulation. 


The numerical results of Ip^crit are shown in Fig. (c). 
When Vp < fp.crit, the effective potential V/ff is maxi¬ 
mized at (sz) = 0, i.e. Kff(l) < 14ff(0). Thus the poten¬ 
tial at the middle point (sz) = 0 becomes a barrier which 
the pseudo spin (sz) cannot pass through [Fig. (a)]. 
This leads to (szit)) > 0 or |/3^(t)| < 1/2 at all the time. 
On the other hand, when the perturbation increases to 
be 1^ > the effective potential is maximized at 

the starting point and > I4ff(0) [Fig. i(b)]. In 

this case, the pseudo spin (sz) can evolve from (sz) = 1 
to (Sz) = —1, and a full transition is reached between 
the two magnetized phases. In general, the maximum 
transition ratio is given by: 


1/3 


2 


max 


A--JA‘^-4V^ 

-’^ben V < Fp,crit, (iQ) 

1, when V > Fp^crit- 


It is clear that |/3niaxP <1/2 when the perturbation is 
below the critical magnitude [Fig. |(d)]. 

The effective theory given by Eqs. ([^ enables a clear 
understanding of the underlying mechanism of the criti¬ 
cal perturbation. The existence of critical external per¬ 
turbation is due to the interaction induced energy mis¬ 
match between magnetized phases and the superposition 
of magnetized phases (stripe state), but not necessar¬ 
ily requiring that the stripe state has an energy higher 
than the magnetized phases. From the effective poten¬ 
tial I4ff((Sz)) we can see that for the regime Eg < Em, 
which implies that ground state of the system is a stripe 
phase, and assuming the system to be initialized in the 
state with a fixed momentum kx = km, the full transi¬ 
tion from the initial phase to another degenerate state at 
kx = —km still necessitates a perturbation exceeding the 
critical value Ip.crit = \Es — i?m|/2. 


Furthermore, the period of the pseudospin oscillation 
{sz{t)) can be solved exactly from the Eqs. (|^. For 
the two different regimes Vp < Ip.crit and Vp > Vp^cnt, 
the oscillation period takes different forms, given by 
(*/"<) 

respectively. Here K{z) is the complete elliptic integral 
of the first kind [SHIj with r< = (A^ — WpY/"^j[2Vp) 
and r> = (41^^ — A^)^/^/A. A more intuitive under¬ 
standing of the pseudospin dynamics can be acquired by 
expanding the oscillation period with the perturbation 
below and over critical value in series of Vp and ^jVp, 
respectively, which gives 

\ i + ^ + W + for E < Ep.erit, 

Tr = { (II) 

TT I -n-A^ I QttA^ I r T/ ^ T/ 

I, Vp 16Vp3 “T 1024 VpS roi V ^ Vp,cnt- 

From the above results we can see that the oscillation 
period increases with Vp when the perturbation is below 
critical value, while it decreases with Vp if the perturba¬ 
tion is over critical value [c.f. numerical results in Fig.j^ 
(d)]. This feature is clearly different from the Rabi os¬ 
cillation in a single-particle two-state system, where the 
oscillation period decreases monotonically as perturba¬ 
tion strength increases. In particular, in the vicinity of 
V « Ip,crit) the period can be approximated by: 


Tr « 5.54 - ^ log(r) -|- - 


[fog(i)+l] 


+ 0{r^), ( 12 ) 


where r = (or r>) and g = A (or 4V),) for Vp > l^.crit 
(or Vp < lp,crit)- The Eq. (12) shows that the oscilla¬ 
tion period diverges logarithmically as the perturbation 
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FIG. 3: Pseudospin spin dynamics governed by the 
effective double-well potential. The initial conditions of 
pseudospin polarization are characterized by (sz(0)) = 1 and 
d{s^(t)) / dt\t^Q = 0. a. The effective double well potential 
maximizes at (sz) = 0 when Vo < Vb.crit, and the pseudospin 
cannot pass through the maximum point (sz) = 0. b. The 
effective double well potential maximizes at (sz) = ±1 for 
Vo > Vb.crit- In this case the pseudospin can pass through 
(Sz) = 0 and thus a full transition occurs, c. The critical 
magnitude of external perturbation as a function of Raman 
coupling Q and interaction energy gsU. d. The pseudospin 
evolution period Tr (red curve) and maximum transition ra¬ 
tio |/3max| (blue curve) versus Raman coupling Q, with the 
interaction energy taken as gsU = O.Sfcp. 


approaches to the critical value [Fig. |^(d)]. On the con¬ 
trary, the oscillation become faster when Vp tunes far 
away from fp.crit- In particular, when 1^ ^ A, the pe¬ 
riod approaches to a “universal” value Tr —>■ tt/A which 
solely depends on interaction energies of the system. On 
the other hand, in the limit of large perturbation V), 3> A, 
the period decreases monotonically to be Tr « tt/V^, im¬ 
plying that the spin dynamics in this regime are domi¬ 
nated by the external perturbation, and are more similar 
as those in the single-particle systems, consistent with 
the numerical results in the previous section. 


Beyond Two-State Resonant Coupling Regime 

Multiple States Correction.- In the previous sec¬ 
tions we have considered only the direct transition be¬ 
tween two magnetized phases induced by the external 


perturbation Vex- On the other hand, it is easy to know 
that Vex can also drive the couplings between the two 
ground magnetized phases and other states. The leading- 
order corrections include the couplings to the two lower 
subband states (6) with momenta kx = ±3fcm, 

and the two upper subband states {6) with mo¬ 

menta kx = ik-m (corresponding to the Hamiltonian 
Hq). Taking into account these additional couplings 
yields a more complicated nonlinear dynamical equation 
which describes the mixing of six components and can be 
solved numerically (see the Supplementary Materials for 
details [35]).. Note that the four additional states are de¬ 
tuned from the ground phases with the energy difference 
per particle AE « 2/cq, which is two orders greater than 
the critical perturbation Ip.crit in the typical parameter 
regime in the experiment. We expect that the density 
of atoms pumped to these states is very small compared 
with the population in the two ground states. As a re¬ 
sult, in deriving the complete dynamical equation, one 
can neglect the interactions within the states (d) 

and (0), while the interactions between them and 
the ground states are still considered [55] . 

The numerical results of the corrections to the maxi¬ 
mum transition ratio and oscillation period, and the cor¬ 
rected critical perturbation are shown with blue curves 
in Fig. 1^ (a-c). We find that, in all the experimentally 
relevant parameter regimes, the inclusion of couplings 
between magnetized phases and other states has negli¬ 
gible corrections to the spin dynamics predicted in pre¬ 
vious sections. In particular, when Vb is tuned away 
from Vb,criti which is the corrected critical perturba¬ 
tion and is slightly smaller than Vb,crit [Fig. (c)], the 
maximum transition ratio |/3maxP is slightly enhanced 
for Vo < I^,crit compared with the result |/3maxP ob¬ 
tained with two-state approximation [Fig. (a)]. Ac¬ 
cordingly, the corrected oscillation period of the spin dy¬ 
namics slightly increases for Vq < Vb.crit and decreases 
for Vb > V^o, crit [Fig. 0 (b)]. These features suggests 
that the leading-order correction strengthens the effec¬ 
tive coupling between the two magnetized phases. This 
phenomenon can be understood that mixing the ground 
phases with the additional states of high energies con¬ 
tributes more intermediate channels to the coupling be¬ 
tween the two magnetized phases, and in consequence 
the net effective coupling (Vg®® = Vb + SVq, with the cor¬ 
rection |dVbl ^ 14) is enhanced. In the experimental 
parameter regime with O.Sfcg and gsU ~ l.Ofcg [32j . 
the inclusion of the excited states alters the oscillation 
period and amplitude (for Vg < 14 . crit) by less than 1%. 
Therefore, we conclude that the two-state approximation 
employed in the previous sections well captures the quan¬ 
tum spin dynamics driven by the external perturbation. 

Momentum Mismatch.- Another approximation 
employed in the previous studies is that we take the 
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Critical 

Perturbation(/(^) 



FIG. 4: Corrections due to the couplings to multiple excited states (blue curves) and momentum mismatch 
(red and green curves), a. Correction to the maximum transition ratio |/3max| ~ |/3max|- b- Correction to the oscillation 
period Tr — Tr. c. The corrected critical perturbation Vo.crit as a function of the Raman coupling strength. The correction for 
Vo.crit with fc_R = fcmin is nearly zero and not observable. The black curve represents the critical value Vb.crit with the two-state 
approximation. Other parameters are taken as gsU = l.Ofco and fl = O.Sfco in the numerical calculation. 


wave vector of standing wave light used to generate 
perturbation to be ku = km- In the experiment, it 
might be difficult to precisely match the wave vector 
kfi with momentum of the magnetized phases, while 
kmin and fco can be precisely determined. The sim¬ 
plest way is to put that kji = kmin or ka = fcg in 
the real experiment, and then the perturbation reads 
14x = lo / cos^ fcija;('!/'|V't + t > 0. From 

Fig. 0 (a,b) we notice that in the typical parameter 
regime {gsU ~ l.Ofcg and 0.2fcg < 0 < 0.6/c§) for mag¬ 
netized phases, the momentum mismatch 6k = k^ — km 
{kji = fcniin or kfi = ko) is a small number, thus we ex¬ 
pect that it only brings minor effects on the quantum 
spin dynamics. Specihcally, due to the momentum mis¬ 
match, the perturbation couples the initial magnetized 
phase with k^ = km to another state at k^ = —km — 2Sk, 
which has a small detuning from the ground magnetized 
state. This effect can slightly reduce the transition ratio 
to be |/3maxl ^ 1 ill fiio regime Vg > Ib,crit- Moreover, 
since the spin polarization and density profile depend on 
the momentum of the condensate, the momentum mis¬ 
match may modify the interaction energies. 

The corrections due to momentum mismatch are 
shown numerically with red and green curves in Fig. 
(a,b) with different perturbation strengths. We find that 
the momentum mismatch slightly affect the oscillation 
periods and the critical value of the perturbation. Es¬ 
pecially, the correction is tiny when the perturbation is 
far away from the critical value, while it is relatively 
more obvious if Vg ~ Vg^crit- This is because in the lat¬ 
ter regime, it takes relatively long time for the system 
to evolve through the intermediate superposition phase 
(stripe state). Note that in such intermediate process 
of passing through the stripe state, interactions play the 
main role in governing the quantum spin dynamics. As a 
result, the aforementioned modifications to the interac¬ 
tion energies due to the deformed spin polarization and 
density prohle may have a more pronounced effect on the 
pseudospin evolution. Nevertheless, we have confirmed 


that in all the typical parameter regimes which are rele¬ 
vant for the experiment, the correction due to momentum 
mismatch is tiny in a large range of Vg. 


Discussions and Conclusion 

The perturbation 14x(r, t) is a conventional spin- 
independent lattice potential applied at t > 0. Differ¬ 
ently from the spin-flip Raman coupling used to creat 
SO coupling, the periodic perturbation 14x can be gen¬ 
erated by a standing-wave light with relatively large de¬ 
tuning, and applying this perturbation should not in¬ 
duce additional heating effect. For the alkali atoms, a 
spin-independent optical dipole potential can be readily 
generated with a 7r-polarized light. 

In experiment, the transition between the magnetized 
phases (pseudospin dynamics) can be detected by mea¬ 
suring the evolution of spin polarization of the conden¬ 
sate according to the relation = (sz(t)) cos20. 

In the typical parameter regimes with gsU ~ I.Ofcg and 
0.2fcQ < O < O.Gfcg, we can verify that km = kmm — ko, 
and cos 26 < 1.0. Thus the pseudospin polarization along 
0 axis nearly equals the spin polarization of the BEC. 

The present study is not restricted by the finite tem¬ 
perature effect, atom loss, or the existence of a weak 
trapping potential. Actually, the finite temperature ef¬ 
fect and atom loss may cause the change in the conden¬ 
sate density n, which effectively modifies the interaction 
energies characterized by gsU and gan. As a result, these 
effects can quantitatively, but not qualitatively affect the 
critical magnitude of the perturbation. Similarly, the 
trapping potential in the real experimental studies may 
also quantitatively affect the predicted threshold of ex¬ 
ternal perturbation. In the presence of a weak trapping 
potential along x direction, the magnetized phase is not a 
state with fixed momentum, but has a small momentum 
distribution around k^ = 6sikm- In this way, the effec¬ 
tive coupling between such two phases can be generically 























modified compared with the case without trapping po¬ 
tential. Nevertheless, the main results predicted in the 
present work will not be changed. 

In conclusion, we uncovered nontrivial quantum spin 
dynamics in a SO coupled Bose-Einstein condensate 
which has two degenerate many-particle ground states 
with opposite finite momenta and being resonantly cou¬ 
pled to each other by an external perturbation. Due to 
the particle-particle interactions, a novel quantum dy¬ 
namical phase transition is predicted that the transition 
of the condensate from one ground phase to the other is 
forbidden when the strength of external perturbation is 
below a critical value, and occurs only when the pertur¬ 
bation exceeds such critical strength. This phenomenon 
is in sharp contrast to the quantum dynamics in a single¬ 
particle two-level system. We developed an effective the¬ 
ory and showed the exact solutions to the quantum spin 
dynamics, which enable a full understanding of the un¬ 
derlying mechanism of the predicted dynamical phase 
transition. It is noteworthy that the present study can 
be generalized to the large spin system with ID synthetic 


SO coupling [3DH12]j in which case the SU{N) or SO{N) 
symmetric interactions of atoms may generically have 
nontrivial effects on the quantum spin dynamics. This 
work opens an important avenue in the study of novel 
quantum spin dynamics in the SO coupled systems with 
interaction, which on one hand can broaden the basic un¬ 
derstanding of new many-body physics and, on the other 
hand, may have interesting applications to the quantum 
control and quantum information science. The high fea¬ 
sibility of the present study will attract experimental ef¬ 
forts in the near future. 
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Supplementary Information 

Appendix A: Derivation of simplified equation of motion 


Pseudo-spin 


The pseudo-spin (s(t)) of the system is defined by (Pauli matrices acting on the pseudospin spanned by the two 
magnetized states): 

((*,), (%), (s.)) = (a- n ((; J) , (“ -■) , (J \)) 

= (23?(ar),23(a/3*),|ap-|/3|2). 

The evolution (i(t)) describes the dynamics of the system. 


Commutation relations 


With the two-state approximation, the evolution of the system can be described by and 

the dynamics of the condensate can be derived according to the equation + Hint + Kx]- The 

commutation relations with respect to the single-particle Hamitlonian and external perturbation are calculated by: 


bPit),Ho] = / (fr 


(SI) 


with = ^ - 




-f ^ ), and 


[V'(t),I4x] = / d^rVbcos^fcox ['0(t),V’’''(r)'0(r)] 


K) 


a(t)(cos6»V’-fe,„t “ sin6>?/>_fc;) 
/3(t)(sin6»V’-fc„t “ cos9il}-ki) 


— + sin 9 cos 9{ai)L + Pi’R)] ■ 


(S2) 
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In Eq. (SI) we have denoted the single-particle Hamiltonian as ho = k^/2 — koUxiJz + i^/2)ax- The commutation 
relations with respect to the interacting Hamiltonian are given by 


bpR^Hint] = ^ d^r 




= y / [cos6»n^(r)?/)^(r) - sm6»n;(r)V’;(r)]-t y / [cos6»n;(r)V't(i’) - sin6»n^(r)V';(r)] 


[Gi + - |/3|2)] - 2Es [3fJ(a/3*) + i^ia/3*)] 


(S3) 


and 


[V'L-dfint] = ^ / d^r + V'l(i')V’l(i')V’t(i')V't(r) 


y / d^r [^V'fl,V’j(r)V’I(r)V't(r)V’TW 


= -^ / d^re [sindn^(r)^/;^(r) — cosdn^(r)^/^^(r)]-|-y / d^re [sindn^(r)^^(r) — cosdn^(r)V’^(r)] 


[Gi - E^(|ap - l/SniV'L - 2Es[n{aP*) - z3(a/3*)]V'i?- 


(S4) 


where Gi = q( 3 s -t ffa), G 2 = j{gs — 5a), 4?^ = G 2 cos^ 2d, and Eg = 2Gi cos^ dsin^ 9. In the last lines of the above 
equations we have applied the mean-field approximation. Then we can now write down the equation of motion in the 
space of tpH and 


. d / a{t) 


= Tdefi 


dt \m 

where the effective two-state nonlinear Hamiltonian is given by 

K) 


I m J ’ 


TdeS — E]^in + y + ^1 + Em,{\a\‘^ — |/3p)s2 -b 2Es [lft(a/3*)sa; — 3(a,5*)sy] -I- VpSx- 


(S5) 


(S6) 


Equation of motion for spin dynamics 


Making use of the substitution a and /3 to (s), the equation of motion becomes 


d(5z} 

dt 


da* 

dt 


■a + c.c. 


Vp{Sy) 


and 

Asy) d{sx) ^ 2 , Ala 

dt dt \dt^ dt 

= 2 ( 32 ) (—A(sj,) + i {A{sx) + Vp )), 

where A = E^ — Em- Therefore, the original dynamical equations can be recast into a set of real differential equations, 
namely, 

{^Y = -2Vp{~Sy), 

{syY = 2{s,){VpEA{sx)), (S7) 

{SxY = -2A{Sz){Sy), 


where (s^)' = d{si)/dt (i = x,y or z). The equation for (sz) can be simplified by substituting the other two equations 
of {sx,y) into it 


(S.)(Sz)"' + 4A2(S2)3(S2)' - (Sz)"(S 2 )' = 0 

+ 2A"(S.)(0)) (S.) - 2A^(Sz)^ = (S 2 )". 


(S8) 
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If the system starts at Sz(0) = ±1, which corresponds to a magnetized phase, the equation of motion can then be 
written as 


(s,)" = (-4F/ + 2A2) (S,) - 2A2(S,)= 

d 


d{sz) 




(S9) 


Furthermore, the other two components of the pseudo-spin are governed by the following equations: 

(S.)' 


(sy) = - 


2K 


P 


(s ) = — ( - V 1 

^ A 'v2(sA) ■ 


(SIO) 


Appendix B: Analytic solutions to the eqnation of motion 


The full solution for pseudo-spin dynamics 


The dynamical equations of the pseudospin can be solved exactly. If starting with the magnetized phase at = km, 
these solutions are 


{Sxit)) 

(Sy(i)) 

iUt)) 


2VpA 


2Vp 

y/iv/-A^ 




cd - A2, > 


(Sll) 


where cn, dn, sd are Jacobi Elliptic Functions [33] ■ In order to solve the period of the pseudospin oscillation, we 
consider two different situations. For Vp < I^,crit) the oscillation period is determined through (s^(t))' = 0. This 

means that \id{t^Wp — A^, = 0, which gives Tr = 2-^^=== AT (i/r<). Here K{z) is the complete elliptic 

integral of the first kind [33], with r< = (A^ — 4I^^)^/^/(2V^). Accordingly, for Vp > Ip.crit, the oscillation period 
is determined through {sx,soi) = 0, i.e. - A^, = 0. This gives that Tr = AT ii/r>), 

with r> = (4V^^ — A^)^'^^/A. Therefore, the period of the system Tr is 


Tr 


fAT 




2iV-n 




K 




iA 


for I^ < I^ 

,crit; 

for Vp > Vp 

,crit ■ 


(S12) 


The power series expansion for K{z) takes the form: 


TT ( 2n\ 


K{z) = < 




2n! 


S V2"" (nir 


0 ^"- for -I < < I 


/C„(z) for < —1, 


(S13) 


where JCn{z) = (log(4|z|) — 272 ^ -I- 27 „) /|zp"+^ and 7 „ = X]”=i * Direct substitutions of the above series yield the 
expansions in the main text. 
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M 



VolVo,cnt 


FIG. SI: Order parameter in dynamic phase transition. The order parameter versus perturbation strength Vo/Vb,crit, 
with Vb.crit = A/2. The graph is independent of the choice of n and 


Order parameter of the dynamical Phase Transition 


We define the time average of {sz{t)) over an oscillation period as the order parameter of the dynamic phase 
transition: 


M{Vo) 


Tr Jo 


K 


Tr 


{szit))dt 

—4V^ 

P 


for Vo < Vo ,crit 1 


0 , 


for Vo > Vo ,crit• 


(S14) 


The critical behavior of the system can be investigated through the expansion of the order parameter M around the 
transition point Vq = Vo^crit- Then we find that 


M = 


TT 

31og2-log(l-2I/p/A) 


for Vo < Vb,crit, 


which shows that the system undergoes a continuous transition when one tunes the perturbation strength Vq through 
the critical point Vb.crit- 


Appendix C: Beyond two-state approximation 
Multiple states correction 


Let ipi and 'ip 2 be the single-particle eigenstates with momentum kx = km and kx = —km in the upper band, 
respectively; t/js idl’i) be the states with momentum 3km {—ikm) in the lower band. We denote O' as the angle 
describing the spin components of the states ipo and i/ 4 . The commutation relations are given by 


[V'3.4,-ffo] = (^^ - \l^klk ^0 + ^'^ ^3.4- 


(S15) 
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For resonant oscillation, the second-quantized form for the perturbation is 


Vex — Vo j ^ (V'I'-2fc™,(TV'feV + '*/’fc'+2fc„(TV'feV^ ■ 


(S16) 


-^k'-\-2km.<y ^ 

By taking into account the additional states, the commutators concerning Vex are now given by 

[ipR, 14x1 = yV'K + y [sin 26 »V’l + cos 26»V'2 + cos {0 - 0') -ips ]; 

[‘<pL,Vex] = yV'L + y [sin ^OipR - cos 20V'i -b cos {0 - 0') -tpi] ; 

[ipi,Vex] = y V'l + y [sin207/>2 - cos20ijjL + sin {0 - 0') ■03]; 

[ 02 , 14 x 1 = y 02 + y [sin 2001 -b cos20ijjR - sin (0 - 0') 04 ] ; 

[ 03 , 14 x] = y 03 + y [cos (0 - 0') 0K -b sin (0 - 0') 0i] ; 

[ 04 , 14x] = y 04 + y [cos (0 - 0') 0L - sin (0 - 0') 02 ] . 

In the calculation, we have ignored all other states which are of higher energies than 0i to 04. This approximation 
can be justified since we shall show that the multiple states correction is indeed negligible. 

The interacting term 


Ffint = V 


^ j d^rga (0|0[[0;0t) -b y (0|0|0t0t + 0|0|0;0;) . 


(S17) 


The commutation relations are 

[0_R, Hint] =Gi ((ipR + aj3* sin^ 200^ -b aj3* sin 20 cos 2002 -b a*/? sin 20 cos (0 — 0') 03) 
-b G' 2 ( cos^ 20 (jap — |/?p) 0fl -b cos 20sin 20 (jop — |/3p)0i) 

[0L, i7int] =G'i(0l -b a*/3sin^ 20ipR — a*/?sin 20 cos 200i -b a/3* sin 20 cos (0 — 0') 04) 
-b G' 2 (cos^ 20 (|/3p — jap) 0i -b cos20sin20 (jap — l/3p)02) 

[01, i7int] =Gi(0i -b a/3* sin^ 2002 — a/3* sin20cos200i -b a*/3 sin 20 sin (0 — 0') 03) 
-b G 2 ( cos^ 20 (1/3p — jap) 01 -b cos 20 sin20 (jap — \/3\'^)iPr) 

[ 02 ,37int] =Gi( 02 -b a*/3 sin^ 200i -b a*/3sin 20 cos 200ij — a/3* sin 20 sin (0 — 0') 04) 
+ G 2 ( cos^ 20 {\af - \pf) 02 + cos 20 sin 20 (jaj^ - y) 

[03, Hint] =Gi( 03 -b a/3* sin20 cos (0 — 0') 0ij -b a/3* sin 20 sin (0 — 0') 0i) 

-b G 2 (cos 20 cos 20' (jap — |/3p) 03) 

[04, i7int] =Gi( 04 -b a*/3 sin20 cos (0 — 0') 0 l — a*/3sin20sin (0 — 0') 02) 

-b G 2 (cos 20 cos 20' (|/3p — jap) 04) 


In the above calculation we have applied the following mean-field approximation 

0|0^ « ^ [jap cos^ 0 -b l/3p sin^ 0 -b cos0sin0 (a/3*e^*^''“' -b c.c)] 
i/i|y « ^ [jap sin^ 0 -b l/3p cos^ 0 -b cos0sin0 (a/3*e^*^''^ -b c.c)] . 


(S18) 


With the above results, we can write down the effective Hamiltonian in the space spanned by the six states including 
two magnetized phases 


H 


( 1 ) 


eff 


^ Eq + Em{Sz) Vp + EsP^ Em 2 {Sz) 

Vp + Eg^* Eo — Ezn{Sz) —Vp2 — Es2^* 

Em2{^z) Vp2 Eg2p, El Ezn{Sz) 

Vp 2 -b Eg 2 ^* Ezn 2 {Sz) Vp -b 33 s^* 

Vpi -b £'s 4^ 0 143 -b Ego^ 

0 14,4 -b Eg 4 p^* 0 


142 + Es2£, 

Em2{Sz) 

+ Egp, 

El -b Em{Sz) 
0 

-^P3 - EgoC 


Vp4 + Eg^C 0 \ 

0 144 -b Eg 4 ^ 

^P3 + EgoC 0 

0 -Vpo - Eg3^ ’ 

E 2 + EjnsiSz) 0 

0 E 2 — EmsiSz) / 


V 
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where the matrix elements read 


El 

Vp 

Vp3 

Em 

Em3 

Es2 

Es4 


kl 




aP*, E, = ^- dklkl + — + ^ + Gi 


2 

V^o 


112 




^ + \ kPkP + -r + TT + G'lj E 2 — —^ - \/ 9 fc 2 fc 2 


02 Vo 




^ sin 20, Vp 2 = — cos 20, 

^ sin (0 — 0'), Vp 4 = ^ cos (0 — 0'), 

G 2 cos^ 20, = G 2 cos 20 sin 20, 

G 2 cos 20 cos 20', Eg = ^ sin^ 20, 

G G 

^ sin 20 cos 20, Ego = sin 20 sin (0 — 0'). 

^ sin 20 cos (0 — 0'). 


Momentum mismatch 


When the wave vector kn of the laser generating the perturbation lattice potential has discrepancy from km, the 
initial magnetized phase with momentum kx = km will be coupled to the state with momentum kx = —km — 5k and 
5k = 2(fcfl — km)- Let f/'o and iph be the state with momentum km and —km — 5k in the lower band, and 9" be the 
angle describing the spin component of the state carrying momentum —km — 5k. The relevant commutation relations 
are 



[pb,Ho\ = -+ 5fc)2A:2 + 

[ipa, Vex] = y ^ Sin (0 + 0") , 

[ipb, Vex] = y (^V'6 + ^ sin (0 + 0") , 

[ipa, Hint] = Gi {ipa + sin^ (0 + 0") a/3*'!/;^) + G 2 ( (cos^29\a\^ — cos 20 cos 20"|/3p) ipa + sin^ (0 — 0") «/?*■;/'&) ) 

[ipb, Hint] = Gi {ipb + sin^ (0 + 0") a* Pip a) + G 2 ((cos^20|/3p - cos 20 cos 20 "|q;P) pjb 

+ sin2 {e-e")a*PtPa). 

The effective Hamiltonian governing the spin dynamics of the present system reads (in the basis of [xpa V’fc]^) 

( 2 ) Eo + EmaWP - ErnbW Vpa + EgaP* \ 

V Vpa+Ega*P Eos + EmalPl^ - Emb\ap J ’ 

where the Hamiltonian matrix elements are 


y = y - Jklkl + ? + 5 + G 


r^O 
2 


^ {kr+5k)‘^ 112 y 

Eos — -^- \j {kr + 5k)^ + y + y + Gi; 


Vp = ^ sin (0 + 0"), Ema = G2 cos^ 20; 

Emb = G 2 cos 20 cos 20", Eg = y sin (0 + 0") + y sin 0 (0 


0 "). 
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